Study on the Rotor Design Method for a Small Propeller-Type Wind Turbine 
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Small propeller-type wind turbines have a low Reynolds number, limiting the number of usable airfoil materials. 


Thus, their design method is not sufficiently established, and their performance is often low. The ultimate goal of 


this research is to establish high-performance design guidelines and design methods for small propeller-type wind 


turbines. To that end, we designed two rotors: Rotor A, based on the rotor optimum design method from the blade 


element momentum theory, and Rotor B, in which the chord length of the tip is extended and the chord length 


distribution is linearized. We examined performance characteristics and flow fields of the two rotors through wind 


tunnel experiments and numerical analysis. Our results revealed that the maximum output tip speed ratio of Rotor 


B shifted lower than that of Rotor A, but the maximum output coefficient increased by approximately 38.7%. Ro- 


tors A and B experienced a large-scale separation on the hub side, which extended to the mean in Rotor A. This 


difference in separation had an impact on the significant decrease in Rotor A’s output compared to the design val- 


ue and the increase in Rotor B’s output compared to Rotor A. 
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Introduction 


In recent years, energy problems due to global warm- 
ing and the exhaustion of fossil fuels have become se- 
rious; thus, the effective utilization of natural energy 
sources is expected. Wind power generation has a rela- 
tively low power generation cost, thereby making it suit- 
able for commercialization. Propeller-type wind turbines 
are currently used for power generation. These wind tur- 
bines have been upsized and are being introduced not 
only on coastal lands and mountainous areas but also on 
offshore locations. On the other hand, the use of small 
propeller-type wind turbines is expected to become fa- 
miliar in environments such as urban areas and suburbs; 


Nomenclature 
a axial induction coefficient 
a’ circumferential induction coefficient 


i ia 2 2 
A wind-receiving area (m^) = 77; 


however, their distribution is insufficient because of lack 
of power generation owing to startup characteristics and 
lower power generation efficiency, damages due to 
changes in wind direction and speed, and noise pollution 
from turbine rotation. 

The blade element momentum theory [1,2] has been 
widely used for the rotor design of propeller-type wind 
turbines; optimum design [3—6] and performance analysis 
[7,8] are conducted. However, a limited number of re- 
ports detail comparisons of optimal properties and flow 
fields of rotors from turbines designed using blade- 
element-momentum-theory-generated design values. Spe- 
cifically, in a small propeller-type wind turbine, the 
Reynolds number is low, and there are limited usable 


Greek letters 

a attack angle (°) 

p blade angle (°) = r,@ Vo 
A tip speed ratio 


B number of blades 

C chord length (m) 

Cw power coefficient = w/ (p Ay}, /2) 
D rotor diameter (m) 

n rotational speed (min’') 

P static pressure (Pa) 

Q torque (N-m) 

r rotor radius (m) 

Re local Reynolds number = wC/v 
v absolute velocity (m/s) 

w relative velocity (m/s) 

W turbine output (W)=Q@ 


airfoil materials; thus, it is important to examine optimal 
properties and flow fields of rotors in the turbine. In re- 
cent years, reverse taper blades, in which the chord 
length is lengthened toward the tip, in contrast with small 
propeller-type wind turbines designed using the blade 
element momentum theory, have been employed, and 
they have been experimentally shown to be effective at 
improving the turbine performance within low rotational 
speed ranges [9,10]. However, the details on the mechan- 
ism of the performance improvement are unknown. At 
this point, a design method for a small propeller-type 
wind turbine has not been established. 

In this research, our final goal was to establish a high- 
performance design guideline and design method for 
small propeller-type wind turbines. Thus, we designed 
rotors using the rotor optimum design method [11] based 
on the blade element momentum theory while designing 
rotors with extended tip chord lengths and linearized 
chord length distributions. Performance characteristics 
and flow fields of these two types of rotors were ex- 
amined with wind tunnel experiments and numerical 
analyses and were compared with their design values. 


Rotor Design Method 


We designed rotors using the rotor optimum design 
method [11], which is based on the blade element mo- 
mentum theory. The summary is presented below [11]. 

Design specifications are as follows: va is the absolute 
velocity at an infinite distance, r, is the tip radius of the 
rotors, r, is the hub radius of the rotors, and B is the 
number of blades. We assumed the tip speed ratio (A) and 
power coefficient (Cy) as follows: 
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kinematic viscosity coefficient (m’/s) 


v 
P fluid density (kg/m*) 

w rotational angular velocity (rad/s) =2 7m/60 
Subscripts 

1 rotor inlet 

2 rotor outlet 

a axial component 

h hub 

r radial component 

t tip 

u circumferential component 

00 infinite distance 


where p is the fluid density and œ is the rotational an- 
gular velocity. Turbine output (W) and rotational speed (n) 
are confirmed using the equations as follows: 


W=} Crom? 6) 
ae 60AVv,, (4) 
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Because design calculations are performed for each 
radius, the blade is subdivided in the radial direction. The 
dimensionless radius (€) is defined as 


gat 
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where r is the radius. In this study, design calculations 
were performed with č = 0.14, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 
0.8, 0.9, and 0.975. Because calculations cannot be per- 
formed when é = 1, € = 0.975 was used for the calcula- 
tion instead. 

The velocity triangle for the arbitrary radius is shown 
in Figure 1 [11]. The angle between the relative velocity 
(w) of the inflow into the blade element and the chord is 
defined as the attack angle (a), and the angle between the 
rotation surface of the rotor and the chord is defined as 
the blade angle (£). The angle between the rotation sur- 
face and the relative velocity (w) is defined as ø. In addi- 
tion, vortices from blades are of arbitrary radius and form 
a spiral vortex surface. At this time, the slipstream was 
not expanded and was assumed to have the same diame- 
ter. We also assumed a condition in which the energy loss 
is minimized. Betz’ law [12] describes where the vortex 
surface flows without deformation. The locus of the vor- 
tex filament present in the vortex surface is shown in 
Figure 2 [11]. Here, w’, is the convection velocity of the 
vortex filament and is perpendicular to the vortex surface. 
In addition, v’ is the apparent convection velocity in the 
axial direction of the vortex surface. We assume the con- 
vection velocity ratio (¢) to be defined by the equation 
given below (¢ is constant regardless of the radius). In 


this research, we used ¢= 0.3 as the initial value. 


t 


C= (6) 
Vv 


o0 


(1+a’)ro 


Axis of rotation 
Rotation surface 


Fig. 1 Velocity triangle 


Fig. 2 The locus of vortex filament in the rotor wake 


The value for ø is obtained from 


par B=] -£j (7) 


The loss coefficient of the blade tip (F) is obtained 
from the Prandtl equation [12] as follows: 


F= Zeus onl 4 | } (8) 


where g, is the g of the blade tip. We used g when ¢ = 
0.975. 
Next, we select the airfoil and obtain the lift coeffi- 


cient (Cz) of the blade element from airfoil characteristics. 


The local Reynolds number (Re) is obtained from 
wC 4zrFcos¢gsindv,,¢ 
ov vC,B 
where C is the chord length and v is the kinematic 
viscosity coefficient. We set the attack angle (a) of blade 
element to obtain the drag—lift ratio (e) using the local 
Reynolds number (Re) and airfoil characteristics by the 
equation as follows: 
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where Cp is the drag coefficient. The power coeffi- 
cient (Cy) can be written as 


Cy =J -J6 (11) 
where J,’ and J,’ are obtained from 
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Sh 
The obtained Jı and J; as well as the assumed power 
coefficient (Cy) are used to calculate ¢ from 


J, -yJ) -4J,Cy 
2J, 

If the calculated ¢ and the ¢ assumed as the initial val- 
ue do not match, we use the ¢ obtained here and recalcu- 
late with o in Equation (7). If ¢ does not converge even 
after repeating this calculation, we change the assumed 
power coefficient (Cw), airfoil, or attack angle (a). If ¢ 
still does not converge, the design specifications must be 
reviewed. 

If ¢ converges, we obtain axial induction coefficient (a) 
and circumferential induction coefficient (a’) from 
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The relative velocity (w) is obtained using the equa- 
tion as follows: 
1- 
pales 


sing (13) 
Chord length (C) is determined from 
C= 47rFcosøsingv é (20) 
wC,B 
Blade angle (£) is obtained from 
p=$-a (21) 


The obtained chord length distribution (C(r)) and 
blade angle distribution (4(r)) determine the blade shape, 
and an optimum rotor can be established. 


Experimental Apparatus and Methods 


Test Rotors 


Schematics of the two types of test rotors used in this 
experiment are shown in Figure 3, and their specifica- 


tions are listed in Table 1. In addition, the chord length 
and blade angle distributions are illustrated in Figures 4 
and 5. Rotor A was designed using the optimal design 
method[11], the above described blade element momen- 


tum theory. Its design specifications are shown in Table 2. 


We set the absolute velocity of the infinite distance (ve) 
as 8 m/s, the rotor tip radius (r,) as 0.25 m, the rotor hub 
radius (r;,) as 0.07 m, and the number of blades (B) as 3 
and assumed a tip speed ratio (A) of 5 and a power coeffi- 
cient (Cw) of 0.4 for the design. For any radial position, 
MELO031 [13] was used as the airfoil, and the attack an- 
gle (a) was 4°. The chord length of the tip was 0 mm. 
However, giving consideration to strength, in Rotor B, 
the chord length was extended such that the maximum 
blade thickness of the tip in Rotor A was 5.0 mm, and the 
chord length distribution was linearized. In addition, the 
blade angle of both rotors was the same. 


Experimental Apparatus and Methods 


An outline of the wind tunnel experimental apparatus 
is shown in Figure 6. This open-type wind tunnel has air 
outlet dimensions of 800 mm x 800 mm. Downstream 
from the center of the wind tunnel outlet, the 1D (500 
mm) position was designated as the standard position, 
and the wind velocity at this position was set as 8 m/s. 
Wind velocity was measured using a Pitot tube and an 


Table 1 Specifications of test rotors 


Rotor A Rotor B 


Tip radius: 7; 0.25m 

Hub radius: r, 0.07 m 

Number of blades: B 3 

Airfoil MEL031 

Tip chord length 0 mm 33.3 mm 


(a) Rotor A 


(b) Rotor B 


Fig.3 Test rotors 
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Table2 Design specifications of Rotor A 
Rotor A 
Flow velocity: v.. 8 m/s 
Tip radius: r, 0.25m 
Hub radius: r, 0.07 m 
Number of blades: B 3 
Tip speed ratio: 2 5 
Power coefficient: Cy 0.4 
Airfoil MELO31 
Attack angle: a 4° 
Wind tunnel Torque meter Tachometer 
Motor 
nN 
Inverter 
a 800 
@ 500 
“750 
E 


Fig.6 Experimental apparatus 


inclined tube manometer. Flow at 1D, a downstream po- 
sition from the wind tunnel outlet, was confirmed to be 
almost uniform. In addition, test rotors were positioned at 
1.25D (750 mm) downstream from the wind tunnel outlet. 
The load of the rotor was controlled using a motor and an 
invertor. Rotational speed (n) and torque (Q) were meas- 
ured with a magnetoelectric rotation detector and a tor- 
que detector, respectively. From these values, the turbine 
output (W) was obtained. Note that the torque obtained 
was corrected by measuring the torque without a rotor. 


Numerical Analysis Methods and Conditions 


We used the general-purpose thermal fluid analysis 
code ANSYS CFX15.0 for the numerical analyses, and 
conducted three-dimensional incompressible steady flow 
analyses. The governing equations are the conservation 
of mass equation [14] and the conservation of momentum 
equation [14]. The SST (Shear Stress Transport) model 
[15] was adopted as the turbulence model. 25°C air was 
used for the working fluid. 

The total computational region is shown in Figure 7. 
The computational region mainly consisted of an external 
region, a middle region and a rotor region. The external 
region was cylindrical with a diameter 10 times the outer 
diameter D of the rotor. Moreover, the lengths upstream 
and downstream were 10 times and 15 times the length of 
the outer diameter D of the rotor, respectively, and were 
measured from the center of the rotor. The numbers of 
computational grids for Rotors A and B were approx- 
imately 4,130,000 and 4,030,000 elements, respectively; 
the total numbers of grids were 6,020,000 and 5,920,000 
elements, respectively. The boundary conditions were the 
flow velocity (vo = 8 m/s) for the inlet boundary, arbi- 
trary rotational speed for the rotor and the static pressure 
(Po = 0 Pa) for the outlet boundary. About the wall sur- 
face boundary, the outer periphery of the external region 
was assigned as the slip condition and wall surfaces were 
assigned as the nonslip condition. Moreover, the bounda- 
ries between the rotating and static systems were joined 
using the frozen rotor [16] technique. 


Results and Discussion 


Comparison of Performance Characteristics 

The relation between the tip speed ratio (A) and the 
power coefficient (Cy) is shown in Figure 8. First, expe- 
rimental and calculated values for Cw of Rotor B were 
compared. Both Cy values reached their maximum at 4 = 
4.5. At this maximum output tip speed ratio (A = 4.5), the 
experimental value was Cw = 0.335, whereas the calcu- 
lated value was Cy = 0.333; these values were in good 
agreement. However, at low tip speed ratios (A = 2.5) or 
high tip speed ratios (A = 6), differences between these 
two values increased. At the maximum output tip speed 


ratio of Rotor B, the experimental and calculated values 
were in good agreement; thus, our calculated results are 
considered to be valid. 
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Fig. 7 Computational region 
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Fig. 8 Correlation between the tip speed ratio and the power 
coefficient 


Next, let us consider the calculated Cy of Rotor A. 
The maximum value of this Cy (0.240) was previously 
obtained using the designed tip speed ratio (A = 5); how- 
ever, this value was significantly lower than the design 
value, Cy = 0.4. On the other hand, the calculated Cy of 
Rotor B had a maximum value that was lower than the 
design value; yet, it was approximately 38.7% higher 
than that of Rotor A. Yet, the maximum output tip speed 
ratio of Rotor B was shifted slightly lower, as described 
previously. We will examine the significant decrease in 
Rotor A’s output compared to the design value and the 
increase in Rotor B’s output compared to Rotor A. 
Comparison of Flow Fields 

The radial direction distributions of the axial compo- 
nents (vg, and va2) and the circumferential components 
(v,; and v,2) are shown in Figures 9 and 10, respectively, 
for the absolute velocity at the inlet/outlet of Rotors A 
and B based on numerical analysis. In both cases, the tip 
speed ratio was A = 5. Figure 9 shows that when the v,; of 
Rotor A and the design value are compared, at any radial 
position, the v,; of Rotor A was significantly larger than 
the design value. On the hub side, the v,; of Rotor B was 
similar to that of Rotor A, but from the mean to the tip 
side, the value of Rotor B was smaller than that of Rotor 


A, and at around 0.8 < r/r; < 0.9, the values were almost 
the same as the design value. However, on the tip side, 
when r/r, > 0.9, the vq, of Rotor B increased rapidly; the 
same was true for Rotor A. On the other hand, v, de- 
creased compared to v,;. Especially in Rotor B, the dif- 
ference between v,; and vy. was large. Therefore, com- 
pared to Rotor B, the flow in Rotor A was completed, 
meaning that the energy of wind was not sufficiently ex- 
tracted. 

Figure 10 shows that v, in Rotors A and B was almost 
0, as designed, and at the rotor inlet, the flow was with- 
out rotation. However, v2 had positive and negative val- 
ues for both rotors, which were especially large on the 
hub side. When comparing the v,2 of Rotors A and B, the 
values were almost the same on the hub side, but from 
the mean to the tip side, the value of Rotor B was higher 
than that of Rotor A. Differences in the v, and v, distribu- 
tion of Rotors A and B are likely to arise from the fact 
that the chord length was mostly the same on the hub 
side but was different from the mean to the tip side. 
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Fig. 9 Axial component distributions of the absolute velocity 
at rotor inlet and outlet (2 = 5) 
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Fig. 10 Circumferential component distributions of the abso- 


lute velocity at rotor inlet and outlet (A = 5) 


The relative velocity (w)) at the rotor inlets of Rotors 
A and B based on numerical analysis, the local Reynolds 
number (Re), and the radial direction distribution of the 
attack angle (a) are shown in Figures 11, 12, and 13, re- 
spectively. In all cases, the tip speed ratio was A = 5. Ac- 
cording to Figure 11, the w; of Rotors A and B was al- 
most the same and agreed well with the design value. 
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Fig. 11 Relative velocity distributions at rotor inlet (A = 5) 
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Fig. 13 Attack angle distributions at rotor inlet (A = 5) 


Therefore, as shown in Figure 12, Re of Rotor A agreed 
with the design value, increasing from the hub to the 
mean side and then decreasing as it approaches the tip 
side. On the other hand, the Re of Rotor B exhibited a 
similar value to that of Rotor A on the hub side but in- 
creased as it approached the tip side, reaching approx- 
imately 0.9 x 10°. This is because the chord length of 
Rotor B from the mean to the tip side was larger than that 
of Rotor A. Figure 13 shows that the a of Rotor A was 
larger than the design value at all radial positions and 
was especially large on the hub side. Meanwhile, the a of 
Rotor B had a similar value to that of Rotor A on the hub 
side but was smaller than that of Rotor A from the mean 
to the tip side, and at 0.8 < r/r, < 0.9 on the tip side, the a 
value became similar to the design value. 

Figures 14 (a) and (b) show limiting streamline on 
suction surfaces for Rotors A and B when 4 = 5. A large- 
scale separation was observed on the hub side of both 
rotors. Especially in Rotor A, the separation extended to 
the mean side. For the MELO031 used in this research, the 
published lower limit of the Reynolds number characte- 
ristic to the airfoil is 1.0 x 10°. Therefore, values for the 
lift coefficient and drag—lift ratio in the Reynolds number 
ranging below 1.0 x 10° were obtained through linear 
function extrapolation. However, because the local Rey- 
nolds number (Re) was extremely low (approximately 0.4 
x 10°) on the hub side of both rotors, the design value 
and actual value for the lift coefficient and drag—lift ratio 
deviated significantly. Thus, the axial component of the 
absolute velocity became large on the hub side, making 
the attack angle (a) much larger than the design value. 
Thus, a large-scale separation was observed on the hub 
side of both rotors. Especially for Rotor A, Re was low 
from the mean to the tip side with a large attack angle (a); 
thus, the separation extended to the mean. The difference 
in the area of the separation likely had an impact on the 
significant decrease in Rotor A’s output compared to the 
design value and the increase in Rotor B’s output com- 
pared to Rotor A. 
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(a) Rotor A (b) Rotor B 


Fig. 14 Limiting streamline on suction surface (A = 5) 


Conclusions 


We examined performance characteristics and flow 
fields of Rotor A (designed using the rotor optimum de- 
sign method based on the blade element momentum 
theory) and Rotor B (the chord length of the tip was in- 
creased compared to Rotor A, and the chord length dis- 
tribution was linearized) through wind tunnel experi- 
ments and numerical analysis. The findings of our study 
are as follows: 

1. In Rotor A, the maximum output tip speed ratio 
agreed with the design value, but the maximum power 
coefficient decreased significantly in comparison with 
the design value. 

2. In Rotor B, the maximum output tip speed ratio was 
low when compared with Rotor A, but the maximum 
power coefficient increased by approximately 38.7%. 

3. A large-scale separation was observed in Rotors A 
and B on the hub side; however, the separation extended 
to the mean side in Rotor A. This difference in the extent 
of the separation is believed to cause the significant de- 
crease in Rotor A’s output compared to the design value 
and the increase in Rotor B’s output compared to Rotor A. 
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